UNCLASSIFIED 

ad  407  784 1 _ 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STAT  ON.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


HOTICB :  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  soy  purpose 
other  than  In  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  cell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


—  X 


Dl  82-0242 


lo 

o 

\c : 


%OE  I NG 


SCI  ENTI  FIC 
RES  EARCH 
LABOR  AT  OR  I  ES 


l 


On  the  Optimum  Transversal  Contour 
of  a  Body  at  Hypersonic  Speeds 


BOEING  SCI^JTIFIC  RESEARCH  LABORATORIES 
FLIGHT  SCONCES  LABORATORY  TECHNICAL  RBORT  NO.  73 


ON  THE  OPTIMUM  TRANSVERSAL  CONTOUR 
OF  A  BODY  AT  HYPERSONIC  SPEEDS 


ANGELO  MIELE  AND  GARY  R.  SAARIS 


MARCH  1963 


z 


tabu  or  caraarrs 


Su— ary 

1*  Introduction 

2.  Dattmination  of  the  drag 
3o  Fomulation  of  the  problaa 

3,1*  Ntctaaary  condition* 

3 .2.  Synmatry  and  Multiplicity  proportion 
4.  Multipit  and  cyan* trie  aolution* 

4.1,  Bodita  of  claaa  Z 

4.2,  Bodita  of  claaa  II 

4.3,  Bodita  of  claaa  III 

4.4,  Diaeuaaion  of  the  reaulta 

3,  General  aolutiona 


Raftrtnoaa 


1 


ON  THE  OPTIMUM  TRANSVERSAL  CONTOUR 
OF  A  BCD!  AT  HYPERSONIC  SPEEDS^ 

by 

ANGELO  and  GARY  R.  SAARIS***** 


SUMARY 

This  papar  refers  to  a  body  whose  length  and  base  area  are  given  and 
considers  the  problea  of  determining  the  transversal  contour  which  mini- 
sizes  the  total  drag  (sub  of  the  pressure  drag  and  the  friction  drag)  under 
the  following  assumptions t  (a)  the  body  is  slender  in  the  longitudinal 
sense)  (b)  the  longitudinal  contour  is  represented  by  a  power  law;  (c)  the 
distribution  of  pressure  coefficients  obeys  Newton's  impact  law;  and  (d)  the 
distribution  of  skin-friction  coefficients  versus  the  abscissa  is  repre¬ 
sented  by  a  power  law.  The  indirect  methods  of  the  Calculus  of  Variations 
are  employed  and  attention  is  focused  first  on  a  class  of  arcs  composed  of 
n  symmetric  cycles  each  covering  the  angular  interval  Zrt/n  (or  2n  segments 
each  covering  the  angular  interval  ri/n )  and,  then,  on  the  more  general  case 
where  these  symmetry  and  multiplicity  properties  are  not  employed.  It  is 

^This  report  supersedes  the  reports  entitled  "A  Study  of  the  Optimum 
Transversal  Contour  of  a  Body  at  Hypersonic  Speeds  Using  Indirect  Methods, 
Part  I"  and  "A  Study  of  the  Optimum  Transversal  Contour  of  a  Body  at  Hyper¬ 
sonic  Speeds  Using  Indirect  Methods,  Part  II". 

(**  ) 

'Director  of  Astrodynanri.es  and  Flight  Mechanics,  Boeing  Scientific 
Research  laboratories. 

(***) 

'Staff  Associate,  Boeing  Scientific  Research  Laboratories* 
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shows  that  the  extremal  configuration  depends  strongly  on  the  friction  pa¬ 
rameter,  a  parameter  which  is  proportional  to  the  ratio  of  the  cubic  root 
of  the  average  friction  coefficient  to  the  average  thickness  ratio* 

Two  distinct  behaviors  are  possible  depending  on  whether  the  friction 
parameter  is  supercritical  (K^  >  1)  or  subcritical  (K^  <  1).  If  the  fric¬ 
tion  parameter  is  supercritical,  the  extremal  solution  is  of  class  I  (a  circle), 
and  the  friction  drag  accounts  for  more  than  two-thirds  of  the  total  drag* 

If  the  friction  parameter  is  subcritical,  the  extremal  solution  is  of  either 
class  II  (any  combination  of  straight-line  segments  tangent  to  a  basic  circle 
whose  radius  depends  on  the  length,  the  base  area,  the  average  friction  co¬ 
efficient,  and  the  exponents  defining  the  longitudinal  contour  and  the  dis¬ 
tribution  of  skin-friction  coefficients)  or  class  III  (any  combination  of 
circular  arcs  of  radius  equal  to  that  of  the  basic  circle  and  straight- 
line  segments  tangent  to  the  circular  arcs);  furthermore,  the  friction  drag 
is  two- thirds  of  the  total  drag*  Among  the  solutions  of  class  II  and 
class  III,  a  particular  set  is  represented  by  starlike  configurations  and 
a  further  subset  by  regular  polygons. 

While  the  solution  is  unique  in  the  supercritical  case,  it  is  not  unique 
in  the  subcritical  case.  More  specifically,  for  each  given  subcritical 
friction  parameter,  an  infinite  number  of  equal-drag  solutions  exist  which 
minimize  the  drag.  For  the  particular  set  of  starlike  configurations,  each 
solution  corresponds  to  a  different  number  of  cycxes;  however,  the  solution 
becomes  unique  if  one  further  constraint  is  imposed  on  the  problem,  such  as, 
the  inner  radius  for  bodies  of  class  II  and  the  outer  radius  for  bodies  of 


class  III 
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Comparison  of  th«  drag  of  extremal  solutions  with  that  of  bodies  of 
revolution  shows  that  the  relative  difference  is  negligible  for  those  sub- 
critical  values  of  the  friction  parameter  which  approach  one  but  is  ap¬ 
preciable  for  those  values  which  approach  zero*  As  an  example*  the  extremal 
configuration  exhibits  %  less  drag  than  the  body  of  revolution  having  equal 
length  and  base  area  for  *  0*8*  25#  less  drag  for  ■  0.6*  and  57#  less 
drag  for  *  0.4. 
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1.  INTRODUCTION 

la  a  previous  paper  (Ref.  1),  the  transversal  contour  of  a  slender 
body  of  given  length  and  base  area  was  determined  in  such  a  way  that  the 
total  drag  (sum  of  the  pressure  drag  and  the  friction  drag)  is  a  minimum. 
Direct  methods  were  employed,  and  the  analysis  was  confined  to  a  body  whose 
cross  section  is  either  a  regular  polygon  or  is  composed  of  a  basic  circle 
external  to  which  are  superimposed  symmetric  segments  of  a  logarithmic 
spiral.  In  this  paper,  these  arbitrary  limitations  on  the  transversal  con¬ 
tour  are  removed,  and  the  minimum  drag  problem  is  investigated  with  the  in¬ 
direct  methods  of  the  Calculus  of  Variations.  The  following  basic  hypotheses 
are  employed:  (a)  the  body  is  slender  in  the  longitudinal  sense;  (b)  its 
longitudinal  contour  is  represented  by  a  power  law;  (c)  the  distribution  of 
pressure  coefficients  is  Newtonian;  and  (d)  the  distribution  of  skin  fric¬ 
tion  coefficients  versus  the  abscissa  is  represented  by  a  power  law. 
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2.  DETERMINATION  OF  THE  DRAG 

Consider  two  aye tens  of  coordinates  (fig*  1):  a  Cartesian  coordinate 
system  x,  y,  z  and  a  cylindrical  coordinate  system  r,  9,  t.  With  regard 
to  the  Cartesian  coordinate  system,  the  z-axis  is  identical  with  the  un¬ 
disturbed  flow  direction,  and  the  xy-plane  is  perpendicular  to  the  z-axiso 
For  the  cylindrical  coordinate  system,  r  is  the  distance  of  any  point  from 
the  z-axis,  and  9  measures  the  angular  position  of  this  point  with  respect 

to  the  zx -plane.  Next,  denote  by  u  ,  u.,  u  the  unit  vectors  of  the  cy- 

r  o  2 

lindrlcal  coordinate  system,  and  observe  that  the  direction  of  u  is  the 
same  as  the  undisturbed  flow  direction.  Also,  denote  by  n  the  unit  vector 
normal  to  the  Infinitesimal  element  of  wetted  area  dS^,  positively  oriented 
outward,  and  by  t  the  unit  vector  which  is  tangent  to  d and  is  in  the 
direction  of  the  local  flow  after  impact.  Consequently,  the  aerodynamic 
drag  per  unit  dynamic  pressure  is  given  by 


where 


?(n  •  u  )2 
2 


(2) 


is  the  pressure  coefficient  associated  with  Newton's  impact  law  and  is 
the  friction  coefficient.  After  it  is  assumed  that  the  tangent  vector  is 


6 


contained  in  the  plane  of  the  vectors  n,  and  it  is  observed  that 


*  *  \  -  Is  x  sl 


dsw 


de  dz 


(3) 


the  aerodynamic  drag  can  be  rewritten  in  the  form 


where  l  is  the  length  of  the  body.  If  the  geometry  of  the  body  is  described 
by  the  equation 


f(rt  0,  z)  ■  0 


(5) 


the  unit  normal  vector  and  the  gradient  of  the  function  f  satisfy  the  re¬ 
lationship 


-  _ZL  fr°r  *  (Vr)  “8  *  V*. 

w '  k  *  <v>2  *  <r 


(6) 


which  implies  that 
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a  •  u 


k  *  <vr>2  *  <r 


n  X  u  ■ 
z 


+  (Vr)< 

fr  +  (Vr)2+ 


ni/2 


n  •  u_  ■ 


r  k  *  (vr>2  *  <r 


(7) 


Consequently,  Eqs.  (4)  and  (7)  lead  to  the  following  relationship  between 
the  aerodynamic  drag  per  unit  dynamic  pressure  and  the  geometry  of  the  body: 


nf 
.  * 


2f> 
z 


{r  *  (Vr>2  *  A 


de  dz  (8) 


which,  because  of  the  slender  body  approximation  (f z/fr)  «  1,  reduces  to 


nf  2f^  i - 


'0  •'0 


215 


(f/r 


d8  dz 


(9) 
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Because  of  hypothesis  (b),  the  longitudinal  contour  of  the  body  is 
represented  by  the  power  law 


f(r,  S,  i)  ir-  (^j  R(0)  ■  0  (10) 

where  m  is  an  arbitrarily  prescribed  exponent  and  R( 0)  is  a  function  which 
describes  the  geometry  of  the  base  contour  (Fig.  1).  Furthermore,  because 
of  hypothesis  (d),  the  distribution  of  friction  coefficients  is  represented 
by  (Ref.  2) 


c,  .  c„(l  .  .)  (ij  01) 

where  Cfa  is  the  average  value  of  the  friction  coefficient  over  the  entire 

length  of  the  body  and  or  is  a  constant  whose  typical  values  are  0  for  the 

idealized  model  in  which  the  friction  coefficient  is  constant,  1/5  for  the 

turbulent  flow  model,  and  1/2  for  the  laminar  flow  model.  In  the  light  of 

(*  ) 

these  hypotheses,  the  aerodynamic  drag  per  unit  dynamic  pressure  becomes 


'  ^From  Eq.  (12),  it  appears  that  the  exponents  defining  the  longitudinal 
contour  and  the  distribution  of  akin-friction  coefficients  must  satisfy  the 
inequalities  m  >  1/2  and  or  <  1. 
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D 

q 


Cfa(1  *  a)i 
■  ♦  l  -  a 


(12) 


where  the  dot  sign  denotes  a  derivative  with  respect  to  the  argument  0* 
The  associated  base  area  is  given  by 


S 


(13) 
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3.  FORMPLATIOH  OF  THE  PROBLEM 


Prior  to  the  fornulatioa  of  the  minima  drag  problem,  it  is  convenient 
to  introduce  the  dimensionless  haute  radius 


P  = 


2(m  +  1  -  Of) 
(2a  -1)0.-  a) 


(14) 


as  well  as  the  drag  parameter  and  the  friction  parameter 


*D 


D  m 

2^?  (2m  -  1)1/3 


2(m  +  1  -  a) 
«  "  *>  Cfa 


(15) 


K* 


(it* 

"(2m  -  1)  (1  -  a)  Cfa' 

w 

2(m  +  1  -  a) 

1 1/3 


the  second  of  which  is  proportional  to  the  ratio  of  the  cubic  root  of  the 
average  friction  coefficient  to  the  average  thickness  ratio*  After  the  re¬ 
lationships  (12)  and  (13)  are  rewritten  as^  ^ 


(*) 

The  subscript  i  denotes  the  initial  point  and  the  subscript  f,  tbs 
final  point. 
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and  after  the  end  conditions  are  expressed  in  the  fora 
0^  ■  0  ,  0j  »  2rr  ,  *  Pj 


(17) 


the  problem  of  minimizing  the  drag  for  a  given  length  and  base  area  is 
formulated  as  follows:  "In  the  class  of  functions  o(  6)  which  are  consistent 
with  the  isoperimetric  constraint  (16-2)  and  the  end  conditions  (17).  find 
that  special  function  such  that  the  integral  (16-1)  is  minimized1' » 

3.1»  Necessary  Conditions 

If  an  undetermined,  constant  Lagrange  multiplier  is  introduced  and 
the  fundamental  function  is  written  as  (Ref.  3) 

F(p,  p,  X)  «  -g  -•  -g  +  2“\/p2  +  f  ♦  \p2 
0  +  P 


(18) 
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it  can  be  recognized  that  the  previous  problem  ie  equivalent  to  that  of 
minimizing  the  integral 


I 


F(p,  p,  X)  d6 

\ 


(19) 


subject  to  the  isoperimetric  constraint  (16-2)  and  the  end  conditions  (17) • 
Hence,  the  Euler-Lagrange  equation  is  given  by 


d 

d0 


S 

ap 


*  o 


(20) 


where 


u  2P 

Vp2*  p2 


(a) 


Since  the  fundamental  function  is  formally  independent  of  the  argument  0, 
this  Euler-Lagrange  equation  admits  the  first  integral 

F  -  #  »  C  (22) 


where  C  is  a  constant  and  where  the  function  f  is  defined  as 


33 


#(p,  p) 


(23) 


This  first  integral  must  be  solved  for  boundary  conditions  consistent  with 
the  transversal! ty  condition 


which  must  be  satisfied  identically  for  any  system  of  differentials  consistent 
with  the  end  conditions  (17).  Therefore,  it  implies  that 


(25) 


Incidentally,  the  solution  may  or  may  not  involve  qpgular  discontinuities 
depending  on  the  values  of  the  average  friction  coefficient,  the  base  area, 
the  length,  and  the  exponents  defining  the  longitudinal  contour  and  the  dis¬ 
tribution  of  skin-friction  coefficients.  Should  angular  discontinuities  be 

(* ) 

present,  the  Erdmann-Weierstrass  corner  conditions 


”  ^The  symbol  &(•••)  denotes  the  difference  between  quantities  evaluated 
after  the  corner  and  before  the  corner. 


T 


(26) 


oust  be  satisfied  at  every  corner  point.  In  particular,  they  iaply  that 
the  integration  constant  C  has  the  same  value  for  every  subarc  composing 
the  extremal  arc.  Finally,  the  application  of  the  Legendre  test  Indicates 
that  the  following  necessary  condition  must  be  satisfied  everywhere  along 
the  arc  Minimizing  the  drag: 


a2F  2o2 
dp  (p  ♦  p  ) 


l  ♦  p 


-•2  2 

3o  -  p 

;  r 1  .rys 

(p  ♦  p  > 


>  0 


(27) 


3.2.  Symmetry  and  Multiplicity  Properties 

Inspection  of  the  functions  F  and  t  appearing  in  the  first  integral 
(22)  shows  that  they  involve  only  even  powers  of  p  and  do  not  involve  the 
argument  0  explicitly.  This  means  that,  if  the  function  p  ■  p(0)  is  a  par¬ 
ticular  solution  of  the  Euler-Ia  grange  equation,  then  the  functions 

p  ■  p(-  9)  ,  p  »  p(0  +  Const)  (28) 


r 

i 

1 

I 

I 

I 

I 

I 

I 


are  also  particular  solutions.  Eq.  (28-1)  expresses  a  symmetry  property, 
since  it  allows  one  to  obtain  from  a  given  extremal  shape  an  infinite  number 
of  equal-drag  shapes  by  drawing  tM  mirror  image  of  the  original  shape  with 
respect  to  any  arbitrary  refer  -  direction.  Eq.  (28-2)  expresses  a 


i 
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multiplicity  property,  sine*  it  allows  one  to  obtain  from  a  given  extreaal 
shape  an  infinite  number  of  equal  drag  shapes  by  rotating  the  original  shape 
through  an  arbitrary  angle  around  the  z-^xls.  Furthermore,  combined  use 
of  Eqs*  (28-1)  and  (28-2)  shows  that  there  exists  a  class  of  variational 
solutions  composed  of  n  identical  cycles  each  of  which  covers  the  angular 
interval  2rt/n  and  is  symmetric  with  respect  to  the  radial  line  joining  the 
origin  of  the  polar  coordinate  system  with  the  peak  point  of  the  cycle • 

Since  these  special  solutions  are  particularly  amenable  to  an  analytical 
approach,  their  properties  shall  be  investigated  in  Section  4.  Then,  the 
more  general  case  is  considered  in  Section 
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4.  MPLTIPI£  AND  SYWETRIC  SOUTTIOMS 

We  shall  now  rastrict  the  analysis  to  ths  class  of  bodies  composed  of 
n  identical  eyeless  by  hypothesis,  each  cycle  involves  a  pair  of  symmetric 
segments,  and  each  segment  does  not  involve  corners.  The  total  number  of 
segments  is  2n;  each  covers  the  angular  interval  rr/n  and  yields  the  same 
contribution  to  the  drag  and  the  enclosed  area.  Consequently,  the  integral 
to  be  extreaized  and  the  isoperimetric  constraint  can  be  rewritten  in  the 
form 


where 


91  ■  0  ,  9f  -  n/n 


(30) 


and  where  the  terminal  radii  are  unconstrained.  Clearly,  the  minimum  drag 
problem  consists  of  minimizing  the  integral 


I  « 


2nF(p, 


p,  X)  d9 


(31) 


subject  to  the  isoperimetric  constraint  (29-2)  and  the  end  conditions  (30) 
with  the  following  understanding:  the  object  of  the  minimization  procedure 
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is  to  find  not  only  the  equation  of  each  segment  but  also  the  number  of 
segments  2n  and,  hence,  the  interval  of  integration  A0  «  rt/n. 

Concerning  the  optimization  of  the  shape  of  each  segment,  the  Euler- 
la grange  equation  (20)  is  still  valid  and  admits  the  first  integral  (22). 
This  first  integral  must  be  solved  for  boundary  conditions  consistent  with 
the  transversality  condition  (24)  which  must  be  satisfied  identically  for 
every  system  of  differentials  consistent  with  the  end  conditions  (30) . 
Therefore,  it  implies  that 


(32) 


and  that 


F, 


(33) 


Furthermore,  if  the  first  integral  (22)  is  integrated  over  the  interval  A0 
and  Eq.  (33)  is  accounted  for,  one  obtains  the  expression 


(34) 


With  regard  to  the  optimization  of  the  number  of  segments,  the  theorem 
of  differentiation  under  the  integral  sign  is  applied,  that  is,  the  derive- 
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tive  of  the  Integral  (31)  with  respect  to  the  parameter  n  is  calculated  and 
the  result  is  set  equal  to  zero*  This  operation  yields  the  relationship 


Sr, 

n  t 


m  0 


(35) 


which,  if  combined  with  Eq.  (3^),  shows  that  the  following  integral  con¬ 
dition  must  be  satisfied  by  the  optimum  contour: 


«d9  -  0 


(36) 


where,  owing  to  the  boundary  conditions  of  the  problem, 

*i  ’  #f  "  0 

From  the  last  two  equations,  it  appears  that  some  fundamental  information 

on  the  nature  of  the  optimum  shape  can  be  obtained  by  investigating  the 

properties  of  the  function  9  in  the  pp-domain  (Fig.  2)*  Clearly,  t  *  0 

(•  ) 

along  the  arc  composed  of  the  euharc 

p  -  0  (38) 


FT 


The  Legendre  condition  (27)  shows  that  Eq*  (38)  is  valid  for  p  s  1  only 
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(• ) 

ending  at  p  >  1  and  the  subarc 

♦(p,  p)  f  p2  +  p2  "  p1*  *  0  (39) 

beginning  at  p  =  1.  Furthermore ,  this  function  is  positive  to  the  left 
of  the  line  *  *  0  and  negative  to  the  right* 

Assuming,  then,  that  the  initial  point  and  the  final  point  are  arbi¬ 
trarily  located  on  the  line  f  ■  0,  three  classes  of  bodies  must  be  analyzed 
(Fig.  2):  (I)  bodies  with  both  end  points  on  the  line  p  *  0,  (II)  bodies 

with  both  end  points  on  the  line  <•  =  0,  and  (III)  bodies  with  the  Initial 
point  on  the  line  p  *  0  and  the  final  point  on  the  line  t«0.  Mathemati¬ 
cally  speaking,  these  classes  are  indicated  as  follows: 


Class 

I: 

^  ■  o  » 

D  • 

H 

o 

Class 

II: 

*i  *  o  , 

♦ 

II 

o 

(**0) 

Class 

III: 

Pi  ■  0  , 

♦ 

►% 

N 

O 

For  each  class,  the  following  subcases  need  to  be  examined:  (a)  the  func¬ 
tion  t  does  not  change  sign  along  the  integration  contour;  (b)  the  function 
f  does  change  sign  and,  therefore,  vanishes  not  only  at  the  end  points  but 
also  at  some  intermediate  points;  and  (c)  the  function  4  is  zero  at  every 
point  between  the  end  points.  The  following  lemmas  are  of  assistance  in 
finding  the  solution: 

The  Legendre  condition  (2?)  shows  that  Eq.  (39)  is  valid  for  p  2  1  only. 
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Lemma  1,  The  function  t  cannot  have  a  constant,  positive  or  negative , 
sign  along  the  integration  contour.  Should  •  have  a  constant  sign,  the 
optimum  condition  (36)  would  be  violated. 

Lemma  2,  The  integration  contour  may  cross  the  line  »  0  at  either 
at  one  point  or  at  infinite  points  but  not  at  any  intermediate  number  of 
points.  Should  crossings  occur,  the  first  integral  (22)  and  Eq.  (39) 
would  be  simultaneously  valid  at  all  points  of  intersection  with  the  impli¬ 
cation  that,  at  these  points. 


Now,  if  the  numerator  and  the  denominator  of  Eq,  (4l)  do  not  vanish  simul¬ 
taneously,  Eq.  (41)  can  be  satisfied  for  only  one  value  of  p.  On  the  other 
hand,  if  the  numerator  and  the  denominator  vanirh  simultaneously,  Eq.  (4l) 
is  satisfied  for  any  value  of  p. 

Applying  the  previous  lemmas,  one  can  find  the  extremal  solution  by 
the  following  exclusion  process]  for  bodies  of  class  I,  contour  (a)  cannot 
be  extremal  because  it  violates  the  first  lemma;  for  bodies  of  class  II, 
contour  (a)  cannot  be  extremal  since  it  violates  both  lemmas,  and  contour 
(b)  cannot  be  extremal  because  it  violates  the  second  lemma;  for  bodies  of 
class  III,  contour  (a)  cannot  be  extremal  since  it  violates  the  first  lemma, 
and  contour  (b)  cannot  be  extremal  because  it  violates  the  second  lemma. 

In  conclusion,  the  extremal  arc  must  necessarily  be  of  type  (c)  and,  there¬ 
fore,  its  most  general  differential  equation  is 
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*  -  0  (42) 

The  particular  differential  equations  associated  with  the  classes  of 
bodies  identified  previously  are  as  followsi 

Class  I:  p  ■  0 

Class  III  ♦  =  0  (43) 

Class  IIIi  p  ■  0  -*-  ♦  =  0 

Regardless  of  the  particular  case,  the  drag  parameter,  the  friction  pa¬ 
rameter,  the  integration  constant,  and  the  Lagrange  multiplier  satisfy  the 
relationship 


(44) 


which  is  obtained  by  combining  Eqs,  (18),  (29),  (34),  and  (36), 

4,1,  Bodies  of  Class  I 

The  contour  of  these  bodies  is  governed  by  the  differential  equation 
(43-1)  which  admits  the  general  integral 

P  «  Cx  (45) 

where  C^  is  a  constant.  Use  of  Eqs.  (29)  supplies  the  parametric  relation¬ 


ships 
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Kjj  -  p(p3  +  2) 


K 


f 


1 

P 


(*♦6) 


which,  after  the  parameter  p  is  eliminated,  lead  to  the  following  minimum 
value  for  the  drag  parameter  (Tig.  3): 


“d 


(47) 


The  associated  integration  constant  and  Lagrange  multiplier  are  given  by 


C 


* 


(48) 


From  the  Legendre  condition,  one  can  see  that  Eqs.  (43)  through  (48)  are 
valid  for  p  <  1  only,  tbat  is,  for  Kj,  s  1.  Thus,  if  one  defines  Kf  ■  1  as 
the  critical  value  of  the  friction  parameter,  the  following  conclusion  is 
reached!  the  extremal  arc  is  a  circle  for  supercritical  values  of  the 
friction  parameter. 


4.2.  Bodies  of  Class  II 

The  contour  of  these  bodies  is  governed  by  the  differential  equation 
(43-2)  whose  general  integral  is 
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p  "  cos(0  +  C2) 


(49) 


where  the  constant  C2*  to  be  evaluated  froa  the  initial  conditions,  is 
given  by 


1 

Pi  "  OOS  C2 


(50) 


EUnination  of  the  constant  froa  these  equations  leads  to  the  result 


P  - - ~jr— -  (51) 

cos  0  -  y  -  1  sin  0 

which  represents  the  straight-line  side  of  a  starlike  configuration  whose 
tendnal  radii  satisfy  the  relationship 


pf  " 


pi 

cos  a  -  Vpi  - 1  5 


(52) 


In  order  to  determine  the  Initial  radius,  the  iso perimetric  condition 
(29 -2)  must  be  combined  with  the  shape  equation  (49)  to  yield 
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(53) 


This  equation,  together  with  Eq.  (50),  supplies  the  relationship  (Fig,  4) 


K 


f 


(54) 


from  which  it  is  clear  that,  for  a  given  friction  parameter,  the  boundary 
value  problem  does  not  admit  a  unique  solution  but  rather  an  infinite  number 
of  solutions,  that  is,  one  initial  radius  for  each  integer  value  of  the 
number  of  cycles*  The  corresponding  final  radius  is  presented  in  Fig.  5 
versus  the  friction  parameter  for  several  values  of  the  number  of  cycles. 

After  the  initial  radius  is  eliminated  from  Eqs.  (51)  and  (54) (  the 
geometry  of  the  optimum  Ghape  can  be  rewritten  in  the  functional  form 

p  ■  p(9,  Kf,  n)  (55) 


which  is  illustrated  in  Figs.  6  and  7*  la  Fig.  6,  the  friction  parameter 
is  held  at  the  constant  value  -  0.778,  and  the  number  of  cycles  is  varied. 
In  Fig.  7,  the  number  of  cycles  is  held  at  the  constant  value  n  *  4,  and 
the  friction  parameter  is  varied. 

It  should  be  noted  that  solutions  of  class  II  exist  if,  and  only  if,  the 
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relationship  (4l)  is  satisfied  everywhere*  Since  this  is  only  possible 
for 


C  -  0  ,  X  -  -  3  (56) 

use  of  Eq.  (44)  yields  the  following  expression  for  the  minimum  drag  pa¬ 
rameter  (Fig.  3): 


h-i  <57 

Since  this  relationship  is  independent  of  the  number  of  cycles  and  the 
initial  radius,  one  concludes  that  all  of  the  extremal  solutions  correspond¬ 
ing  to  a  given  value  of  the  friction  parameter  have  the  same  drag.  Further¬ 
more,  because  of  Eqs.  (29-1)  and  (42-2),  each  solution  is  such  that  the 
friction  drag  is  two- thirds  of  the  total  drag*  Finally,  from  the  Legendre 
condition,  one  can  see  that  Eqs*  (49)  through  (57)  are  valid  for  p  >  1  only* 
If  this  statement  is  combined  with  the  isoperimetric  condition  (29-2)  and 
the  mean  value  theorem  is  applied,  one  concludes  that  star like  solutions 
of  class  II  are  possible  for  <  1  only,  that  is,  for  subcritical  values 
of  the  friction  parameter. 

Polygonal  solutions*  The  present  starlike  solutions  contain  polygonal 
solutions  as  a  particular  case.  For  oi  *  1,  Eq*  (51)  reduces  to 


p  “  coTT 


(58) 
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which  is  the  equation  of  the  side  of  a  regular  polygon  cireuaacribed  on 
the  circle  pal.  The  number  of  sides  (i.e.  t  cycles)  of  the  polygon  is 
the  lowest  among  all  the  starlike  solutions  of  class  II  corresponding  to 
a  given  value  of  the  friction  parameter  and  is  to  be  determined  from  Eq.  (54) 
which  reduces  to  (Fig.  8) 


Kf  -  VS cot  S  (59) 

Clearly,  these  polygonal  solutions  are  physically  realistic  only  for  a  dis¬ 
crete  set  of  values  of  the  friction  parameter,  that  is,  those  values  for 
which  n  is  an  integer.  Incidentally,  the  number  of  sides  increases  mono- 
tonically  with  the  friction  parameter.  In  particular,  for  Kf  1,  the 
number  of  sides  becomes  infinitely  large,  that  is,  the  extremal  solution  is 
a  circle.  On  the  other  hand,  for  -*■  0,  the  number  of  sides  tends  to  two, 
that  is,  the  extremal  solution  is  a  flat  plate. 

4.3.  Bodies  of  Class  III 

The  contour  of  these  bodies  is  governed  by  the  differential  equations 
(43-3)  and  involves  two  subarcs t  one  is  circular  and  the  other  is  a  straight 
line  tangent  to  the  circle.  Because  of  continuity  considerations,  the  transi¬ 
tion  point  from  the  circle  to  the  straight  line  is  located  at  p  «  1,  which 
means  that,  at  this  point,  the  two  subarcs  are  tangent  to  one  another. 
Analytically,  the  extremal  arc  is  represented  by 
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pal  ,  0  <  8  <  C 

p  *  cos(K"y  •  « s  0  s  Vn 


(60) 


where  «  is  the  angular  interval  corresponding  to  the  circular  portion*  In 
particular*  the  final  radius  is  given  by 


pf 


(61) 


In  order  to  deteraine  the  unknown  angle  c,  the  i soper imetric  condition 
(29-2)  must  be  combined  with  the  shape  equation  (60)  to  obtain  (Fig.  9) 


d© 

— 277” 

cos  (0  - 


(62) 


Elimination  of  the  parameter  e  from  the  last  two  equations  yields  the  re¬ 
lationship 


K 


f 


1/2 


n/n 

-  arc  tan 


(6?) 


which  supplies  the  final  radius  in  terms  of  the  friction  parameter  and  the 
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n labor  of  cycles  (Fig.  5).  For  a  given  friction  parameter,  the  boundary 
value  problem  does  not  admit  a  unique  solution  but  rather  a  finite  number 
of  solutions,  The  particular  solution  e  •  0  corresponds  to  a  regular  poly¬ 
gon  and,  as  such,  is  governed  by  Eqs.  (58)  and  (59)*  The  remaining  eolu- 
tions  are  characterized  by  a  number  of  cycles  less  than  that  of  the  regular 
polygon. 

After  the  parameter  c  is  eliminated  from  Eqs,  (60)  and  (62),  the  geometry 
of  the  optimum  shape  can  be  rewritten  in  the  functional  form  (55)  which  is 
illustrated  in  Figs.  10  and  ?,  In  Fig.  10,  the  friction  parameter  is  held 
at  the  constant  value  *  0,95,  and  the  number  of  cycles  is  varied.  In 
Fig.  7,  the  number  of  cycles  is  held  at  the  constant  value  n  «  4,  and  the 
friction  parameter  is  varied.  Incidentally,  the  minimum  drag  parameter  is 
still  represented  by  Eq.  (57)  and  is  plotted  in  Fig.  5.  Since  this  re¬ 
lationship  is  independent  of  the  number  of  cycles  and  the  final  radius,  all 
of  the  extremal  solutions  which  correspond  to  a  given  value  of  the  friction 
parameter  have  the  same  drag.  Furthermore,  each  solution  is  such  that  the 
friction  drag  i6  two-thirds  of  the  total  drag.  Because  of  the  Legendre 
condition,  these  solutions  are  valid  for  p  s  1  only.  If  this  statement  is 
combined  with  the  isoperimetric  condition  (29-2)  and  the  mean  value  theorem 
is  applied,  one  concludes  that  starlike  solutions  of  class  III  are  possible 
for  Kj  1  only,  that  is,  for  subcritical  values  of  the  friction  parameter. 

4.4.  Discussion  of  the  Results 

From  the  previous  analysis,  it  is  evident  that  the  friction  parameter 
has  a  considerable  influence  on  the  nature  of  the  extremal  solution.  If 


29 


one  defines  Kf  »  1  as  the  critical  value  of  the  friction  parameter,  two 
distinct  behaviors  are  possible  depending  on  whether  the  friction  parameter 
is  supercritical  (K^  >  1)  or  subcritical  (K^  <  1).  If  the  friction  pa¬ 
rameter  is  supercritical,  the  solution  is  of  class  I  (circular),  and  the 
friction  drag  contributes  more  than  two-thirds  of  the  total  drag.  If  the 
friction  parameter  is  subcritical,  the  solution  is  of  either  class  II  (star- 
like  with  straight-line  segments)  or  class  III  (starlike  with  circular  arcs 
p  «  1  alternating  with  straight-line  segments  tangent  to  the  circular  arcs). 
The  polygonal  solutions  are  a  subset  of  either  class  II  solutions  or 
class  III  solutions  which  occur  only  for  discrete  values  of  the  friction 
parameter. 

While  the  solution  is  unique  in  the  supercritical  case,  it  is  not  unique 
in  the  subcritical  case.  More  specifically,  for  each  given  value  of  the 
friction  parameter,  there  exists  an  infinite  number  of  extremal  solutions 
of  class  II  and  class  III.  However,  the  solution  becomes  unique  if  some 
further  constraint  is  imposed  on  the  problem,  that  is,  the  initial  radius 
for  bodies  of  class  II  and  the  final  radius  for  bodies  of  class  III. 

The  solutions  of  class  II  as  well  as  those  of  class  III  are  characterised 
by  the  same  value  of  the  minimum  drag;  furthermore,  the  friction  drag  is 
two-thirds  of  the  total  drag.  Hence,  the  drag  of  these  solutions  is  less 
than  that  of  the  equivalent  body  of  revolution  (subscript  R),  that  is,  the 
body  of  revolution  having  equal  length  and  base  area  (Fig.  3)»  Specifically, 
the  drag  ratio  is  given  by 


D 


*4 


(64) 
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and  is  plotted  la  Fig.  (11)  versus  the  friction  parameter.  Clearly,  the 
drag  reduction  associated  with  the  use  of  starlike  cross  sections  is  negli¬ 
gible  for  values  of  near  unity  but  may  be  appreciable  for  values  of 
near  zero.  Its  magnitude  is  5#  for  =  0.8,  2%  for  ■  0.6,  and  57#  for 
Kf  *  0.4. 


31 


5.  GENERAL  SOLUTIONS 

The  next  step  In  the  analysis  is  to  remove  the  limitations  of  the 
previous  section,  that  is,  to  find  the  extremal  arc  without  the  restriction 
that  the  solution  be  necessarily  composed  of  n  identical,  symmetric  cycles* 
Hence,  the  problem  is  to  find  solutions  for  the  first  integral  (22)  sub¬ 
jected  to  the  end  conditions  (25)  and,  in  the  case  where  discontinuities 
occur,  to  find  the  corner  conditions  (26).  While  a  straightforward  analyti¬ 
cal  procedure  to  find  the  new  optimum  has  not  been  found,  it  has  been 
observed  that  the  differential  form  (42)  of  the  solution  for  multiple,  sym¬ 
metric  bodies  is  also  valid  in  the  most  general  case.  The  proof  of  this 
statement  lies  in  the  fact  that  any  shape  governed  by  the  equation  #  •  0 
satisfies  (a)  the  end  conditions  (25),  (b)  the  corner  conditions  (26),  and 
(c)  the  first  integral  (22)  providing  the  fundamental  function  satis'ies 
the  relationship 


F  =  C 


(65) 


everywhere  along  the  extremal  arc.  It  can  be  verified  that  each  of  the 
three  classes  of  shapes  described  by  Eq6.  (4j)  is  consistent  with  Eq.  (65) 
providing  the  integration  constant  and  the  lagrange  multiplier  be  defined 
by  Eqs.  (48)  in  the  supercritical  case  and  by  Eqs.  (56)  in  the  subcritical 
case. 

In  conclusion,  if  the  friction  parameter  is  supercritical,  the  extremal 
solution  is  still  a  circle;  furthermore,  the  minimum  drag  parameter  obeys 


Eq.  (47).  On  the  other  hand,  if  the  friction  parameter  is  subcritical.  any 
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combination  of  circular  arcs  o  ■  1  and  straight  lines  tangent  to  th«  circu¬ 
lar  area  is  aa  extremal  (see  Fig.  12  for  none  typical  examples)';  furthermore, 
the  minimum  drag  parameter  obeys  Eq.  (3?)* 

In  closing,  a  word  of  caution  is  needed  with  regard  to  the  assumed  dis¬ 
tribution  of  pressure  coefficients.  Since  Newton's  law  is  an  empirical 
law  which  was  found  to  be  experimentally  valid  for  convex  bodies  at  moderate 
hypersonic  Mach  numbers,  the  results  of  this  report,  are  satisfactory  only 
for  the  circular  solutions,  the  polygonal  solutions,  and  the  starlike  so¬ 
lutions  of  class  III.  For  the  starlike  solutions  of  class  II,  the  ac¬ 
ceptability  of  the  results  decreases  as  the  concavity  of  the  extremal  in¬ 
creases.  Therefore,  these  solutions  are  indicative  of  trends  only;  in  par¬ 
ticular,  the  associated  aerodynamic  drag  should  be  recalculated  with  more 
sophisticated  fluid  mechanics  techniques  than  those  employed  here. 
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Fig.  1.  Coordinate  systems 
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Fig.  2.  Alternative  integration  paths  in  the  pp-plane. 


Fig.  3.  Minimum  drag  parameter. 


Fig.  4.  Initial  radius  for  bodies  of  class  II. 


Bodies  of  class  H  and  class  IQ,  n  =  4 


Fig.  7.  Examples  of  solutions  of  class  II  and  class  III. 
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Fig.  9.  Angular  interval  corresponding  to  the  circular  portion 
of  bodies  of  class  III. 


Bodies  of  class  IE,  Kf  — 0.93 
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Fig.  12.  Examples  of  solutions  of  class  II  and  class  III. 


